Dynamical regimes of a quantum SWAP gate beyond the Fermi golden rule. 
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We discuss how the bath's memory affects the dynamics of a swap gate. We present an exactly 
solvable model that shows various dynamical transitions when treated beyond the Fermi Golden 
Rule. By moving continuously a single parameter, the unperturbed Rabi frequency, we sweep 
through different analytic properties of the density of states: (I) collapsed resonances that split at an 
exceptional point in (II) two resolved resonances ; (III) out-of-band resonances; (IV) virtual states; 
and (V) pure point spectrum. We associate them with distinctive dynamical regimes: overdamped, 
damped oscillations, environment controlled quantum diffusion, anomalous diffusion and localized 
dynamics respectively. The frequency of the swap gate depends differently on the unperturbed Rabi 
frequency. In region I there is no oscillation at all, while in the regions III and IV the oscillation 
frequency is particularly stable because it is determined by the environment's band width. The 
anomalous diffusion could be used as a signature for the presence of the elusive virtual states. 

PACS numbers: 03.65.Yz, 03.65.Xp, 42.50.Dv 
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I. INTRODUCTION 



It is well known that processing information using 
quantum mechanics makes possible communication pro- 
cedures and computational tasks that could outperform 
classical devices in terms of security or speed Q . There 
are extensive experimental efforts in many fields to re- 
alize implementations of the necessary building blocks 
for quantum information processing. The key ingredient 
is an externally controlled evolution of a superposition 
state which implies a form of quantum parallelism. The 
alternatives range from electrons spins ^ , superconduct- 
ing circuits [1, U, quantum electrodynamic cavities 0, Q , 
and optical ion traps to nuclear magnetic resonance 
(NMR) in liquids % and sohd state 0. One of the most 
important building blocks is the swap gate, where a sys- 
tem S jumps between two degenerate states, A and B, 
when the coupling Vab is turned on. Starting on state A, 
the return probability oscillates with the Rabi frequency 
loq = 2VAB/fi- However, it is clear that it is not feasible 
to isolate completely any real system. In practice, the 
interactions with an environment £ fld,^ perturb the 
evolution, smoothly degrading the quantum interferences 
with a "decoherence" rate, l/r^. This rate is usually 
identified with the S£ interaction rate l/r^g, typically 
evaluated from the Fermi Golden Rule (FGR). In spin 
systems, A and B could be the tJ, and if spin configura- 
tions respectively. Weak interactions (l/r^g ^ 2^0 ) pro- 
duce a slightly slower oscillation which decays at a rate 
1/t0 = l/(2T^g). There are experimental conditions, 
however, where the observed frequency shows a dynam- 
ical transition on its dependence of the S£ interaction 
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[9|. In fact, the swapping frequency is a non-analytic 
function of the interaction rate. At a critical strength 
l/rjg = 2ti)o, the oscillation freezes indicating a transi- 
tion to a new dynamical regime. The initial state now de- 
cays to equilibrium at a slower rate 1 /t^ oc ujq t^^ , which 
vanishes for strong S£ interaction. This last regime can 
be seen as a quantum Zeno phase, where the dynamics is 
inhibited by the frequent "observations" of the en- 
vironment. Such quantum freeze can arise as a pure dy- 
namical process, governed by strictly unitary evolutions 
(see Refs. [l^ EH)- Indeed, some of the phenomenol- 
ogy of that transition was not foreign to spectroscopists. 
The collapse of the independent resonance lines leads to 
the exchange and motional narrowing addressed by Van 
Vleck Ts'] and Bloembergen, Purcell, and Pound in 
the 1940's and synthesized in the analytic properties of a 
phenomenological classical probabilistic model by P. W. 
Anderson Jl7|. The quantum description of the phase 
transition [9| required a self-consistent calculation of the 
oscillation in presence of the S£ interaction. While the 
system remembers its previous state, the environment 
was described in the "fast fluctuations" approximation 
where it retains no memory of its previous state . 

In this work we present a simple and exactly solvable 
quantum model that, while only describing the coherent 
part of a spin SWAP gate dynamics [l^, presents the 
dynamical transitions. Our goal is to deepen the under- 
standing of how bath's memory affects the system dy- 
namics beyond the Fermi golden rule. A complementary 
vision for the dynamics under the action of a given Hamil- 
tonian is, of course, the spectral representation which can 
be studied as function of the S — £ interaction strength. 
However, the energy representation hides much useful dy- 
namical information in subtle spectral properties, such 
as resonances that collapse at the "exceptional points" 
(EPs) in the complex plane. Other unusual properties 



involve resonances that shrink and jump into the non- 
physical Riemann sheet to become virtual states, which 
sometimes are also called antibound states [20, [2i|. Ulti- 
mately, these resonances can transform themselves into 
isolated singularities on the real axis, accounting for lo- 
calized states. In our model all of these transitions will 
appear naturally through the variation of a single con- 
trol parameter. There are a number of physical systems 
that show some of these delicate spectral properties: The 
EP or collapse of resonances has been observed in crys- 
tals of hght HI], electronic circuits [1^, propagation of 
light in dissipative media 24, "25^, vacuum Rabi splitting 
in semiconductors cavities ^26i] . in microwave billiards 
[27! . m, |2^, [lOl , detuned but coupled piano strings [3l| , 
and there are a number of examples drawn from electron- 
paramagnetic resonance [s^l , and solid state NMR ^] . It 
also may appear in many theoretical models: e.g., de- 
scribing the decay of superdeformed nuclei 'ss'l, phase 
transitions and avoided level crossings [34, 35, 36], cou- 
pling of bound states to o pen decay channels [33|, geo- 
magnetic polarity reversal [38| , tunnelin g b etween quan- 
tum dots [3^, |40|. optical microcavity [4ll |. vibrational 
surface modes [42|, and in the context of the crossing 
of two Coulomb blockade resonances [i^ . On the other 
hand, the virtual-localized transition has been addressed 
in the context of the famous n — p singlet system ^2^ , 
models of stabilization of quantum mechanical bin ding 
by potential barriers [i^l, Feshbach resonances lisl . [ig j 
and stability of atomic and molecular states [13, 13 ■ 

This paper is organized as follows. In Sec. II the model 
is presented. It is solved in Sec. Ill using a Green's 
function formalism |49] and its analytical properties are 
analyzed. In Sec. IV the different parametric regions are 
associated with the analytic properties of the density of 
states: region I) collapsed resonances; region II) resolved 
resonances; region III) out of band resonances; region 
IV) virtual states; and region V) pure point spectrum. 
These parametric regions are associated with distinctive 
dynamical regimes: overdamped decay, damped oscilla- 
tions, environment controlled quantum diffusion, anoma- 
lous diffusion and localized dynamics, respectively. The 
anomalous diffusion then could be considered as a signa- 
ture of virtual states. The summarizing conclusions are 
presented in the last section. 
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FIG. 1: Representation of the system and the environment. 
The first two sites are connected through Vab- The hop- 
ping amplitude inside the environment is V. The system- 
environment interaction is Vo- 



where 

Hs = Ea \A) {A\ + Eb \B) {B\ - Vab {\A) {B\ + \B) {A\) 

(2) 
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and with n > 1, is the state localized at the nth site 
of the chain. Vab, Vq and V are the positive hopping 
amplitudes between two contiguous sites. The first two 
sites, linked by Vab, are what we call the system S. The 
rest of the chain is what we take as the environment £. 
In this work, we are interested in the case Ea — Eb — 
En = 0. This choice of parameters together with the 
asymmetric connection to the environment enables a rich 
spectrum. 

This model applies to implementations of NMR quan- 
tum information processors by resorting to the Jordan- 
Wigner transformation |9i| ^5Q] that maps spins into a 
Fermionic system. More specifically, spin chains with 
XY interaction are transformed, at the experimentally 
relevant high temperatures, into a system of indepen- 
dent fermions in a tight binding chain [5l| . In particular, 
Ref. |9i] shows how the Swap Gate in solid state NMR, is 
obtained from a coherent dynamics as that of our model 
plus an incoherent contribution. The fermionic model, 
having incorporated the temperature in the appropriate 
regime, is a natural alternative to the various spin-boson 
models, where a number of transitions are also known 
to exist as a function of temperature and interaction 
strength [s^ ]. 



II. THE MODEL 

Our model, schematized in Fig. [1] is a tight-binding 
semi-infinite linear chain. The Hamiltonian has one part 
that describes the two states of the system, another that 
describes the infinite degrees of freedom of an environ- 
ment and a third part that is the interaction between 
both. It is expressed as 
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III. ANALYTIC SOLUTION 

In this work we use the Green's Function (GF) formal- 
ism to solve the model. Starting with the isolated system 
Hs, we define the GF matrix G, 



(e) = (H5 - £l)" 



(5) 



In order to make a fully consistent definition of Gn,m(£) 
as the Fourier transform of the retarded propagator 
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Gn,m{i) we assume that each site n has an intrinsic decay 
process i?„ — > En — ipn- This solves the GutzwiUer objec- 
tion on this aspect [53|. However, since these imaginary 
parts are considered infinitesimal, they are not written 
explicitly in what follows. 

The first diagonal component of G (e) gives us infor- 
mation about the system dynamics. This component be- 
comes 



and 



G^A is) 



1 



e — 



(6) 



Here, we can see that the poles of G^^\ (e) give the sys- 
tem bonding and antibonding energies Ei^2 = ^Vab- 
We consider an isolated system of two sites and excite 
the first one |1,0) = \A). The system evolves oscillating 
between the states \A) and |0, 1) = \B) with a character- 
istic Rabi frequency u>o = . This oscillation is used 
to generate a SWAP gate by letting the Hamiltonian act 
during a time iswap = 

When we take into account the environment, the first 
diagonal component of the GF becomes 



Gaa (e) 



1 



^AB 



(7) 
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e>2V 

\e\ < 2V 
e < ~2V. 



(12) 



A brief commentary on the complex self-energies, i.e., 
non-Hermitian terms, is in order. Its appearance, either 
in a FGR calculation or in the exact solution of Eq. [9l re- 
lays on the fact that the new eigenstates are completely 
orthogonal to the unperturbed ones. In our case they 
are extended states enabled by the consideration of the 
thermodynamic limit: the number of states becomes in- 
finite before 77 is allowed to reach [59]. As discussed 
above, adding this small imaginary part means putting 
the system in contact with an additional environment £' . 
In this situation, Eq. [9] produces a self-energy containing 
a square root function of the energy instead of a ratio 
among polynomials that results when Eq. \E\ is applied 
to a finite system. The sign in front of the square root 
in Eqs. [TT] and [12] is chosen to ensure the physical (i.e. 
decaying) behavior when 77 < 0. As discussed in Ref. 
[soj the presence of non-Hermitian terms is fundamental 
in allowing a dynamical phase transition. In the general 
context of non-Hermitian quantum mechanics a similar 
conclusion holds [60j . 

The poles of Eq. [71 control the dynamics of the system 
and can be obtained analytically as 



Typically the self-energy E ~ A — iF is evaluated within a 
Fermi golden rule approximation [HJ as an e independent 
complex number. Since this would imply neglecting all 
dynamics and memory effects of the environment , such 
procedure could miss some subtle behaviors jH, [H, . 
Our model enables the evaluation of the exact self-energy 
of an environment represented by the semi-infinite chain, 
and hence accounts precisely for these "memory effects" . 
According to the continued fractions solution ^] of the 
renormalized perturbation expansion [ssj . 



S(e) = 



e-I](£) 
which sums up to the form 

I](e) = A(e)-ir(e), 



(8) 



(9) 



(10) 



(13) 

Equation [T51 has four solutions, see Fig. [3 When they 
have imaginary parts, only those with negative ones rep- 
resent a decaying response to an initial condition. These 
imaginary parts are precisely the exponential decay rate 
in the self-consistent Fermi golden rule [6l|. When the 
four poles are real, the physical ones should approach to 
the isolated system poles (shown in Fig. [2| with dashed 
lines). The real and imaginary parts of the poles are 
shown as a function of the system hopping Vab for fixed 
values of Vb = 0.8^. The solutions not satisfying the 
above conditions are indicated with dotted lines. Here- 
after, we will refer as "the Poles" of the GF only to those 
indicated by the continuous line. 

The Local Density of States (LDOS), can be evaluated 

as 



NA{e) = -- lim lmGAA{{En 

TT r;^0+ 



i»7«}v„-£)- (14) 



with 



A(e) 



(f) -V2 ^>2F 

I kl<2y 



(11) 



This definition is equivalent to the standard one, 
— i lim^^o+ Im (e + i?y) , for most practical purposes. 

An advantageous feature of the present model is that 
the LDOS for lei < 2V can be factorized as 



Na (e) = ^1 {e) X Li (e) x (e) . 



(15) 
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FIG. 2: (Color online) Real (a) and imaginary (b) part of 
the Green function poles vs. Vab for Vb = 0.8. All energies 
and hoppings are in V units, (a) The dashed line represent 
the poles for an isolated system. Different colors identify the 
poles. The nonphysical ones are represented with dotted lines. 
The vertical dotted lines divide the dynamical regions. 



Here, Li and X2, are Lorentzian functions (LFs), 
A^i,is the density of directly connected states (i.e. 
LDOS of the first site of the semi infinite chain). 
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The LFs Li (e) and L2 (e) are related with the real and 
the imaginary part of the GF's poles. Their centers move 
with the real part of the poles, and their widths are de- 
termined by the imaginary parts. 

In Fig. [2] it is observed that there are regions with dif- 
ferent analytical behaviors, some of them separated by 
abrupt changes that are consequence of the nonanalyti- 
cal behavior of the GF poles. The difference between the 
physical poles has a real part uj — |Re (e^i) — Re (£^2)! = 
2 |Re (£ri)|, representing an effective Rabi frequency. The 
imaginary part is associated with the decay rate into 
environment's states. This work mainly focuses on the 
study of the behavior of w, as the parameter that char- 
acterizes the dynamics. 

The relation between GF's analytic properties and the 
dynamics is clarified by writing the survival probability 



FIG. 3: (Color online) Zoom of Fig. [2] in the region where 
virtual states and out of band resonances occur. All energies 
and hoppings are in V units. 



in the energy-time representation 
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PAA{t) = J de J ^GAAie + ^hij) 



— 00 —00 



xG\^{e — \hjj) exp (-icjt) 
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(17) 
(18) 



This function measures the probability to find a particle 
in the site A at time t, provided that the system has had a 
particle at the same site at time t = 0. When the system 
is isolated, Paa {t) oscillates with frequency luq, which 
coincides with uj. When the environment is taken into 
account Paa (t) evolves in a more complex way. In spite 
of this complexity, the evolution at short times can be 
described by an exponentially decaying oscillation with 
frequency uj. In the next section we will present a deeper 
analysis of all spectral regions and their main dynamical 
characteristics. 



IV. PARAMETRIC REGIONS 

By analyzing the qualitative features of the LDOS as 
function of the control parameter Vab we find five para- 
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metric regions separated by well defined critical values. 
We enumerate them from I to V as they are appearing 
by increasing Vab and name them according to the main 
features in the LDOS and the behavior of the GF's poles 
in the complex energy plane. The most common situa- 
tion occurs when the two states of the isolated system 
are mixed with the environment continuous and hence 
acquire a finite mean life. This is region II of resolved 
resonances, a regime typically described by the FGR. 
When the interaction with the environment becomes 
strong enough, at the exceptional point appears the non- 
trivial transition to Region I of collapsed resonances. 
This is the regime where exchange narrowing occurs (17} . 
In the other extreme, we may consider that the internal 
interaction of the system is much stronger than the band- 
width of the environment's continuous spectrum. Hence, 
the system's bonding and anti-bonding states are pushed 
away from the band according to perturbation theory 
and they will remain localized. This is Region V of pure 
point states. Region III with out-of-band resonant 
states and Region IV of virtual states, would go almost 
unnoticed unless the internal and external interactions 
are very similar, see Fig. O In these cases the separation 
between system and environment becomes very delicate. 

The main features of these regions can be identified 
by following the GF's poles into the complex plane as 
represented in Fig. [5] and zoomed in Fig. [31 We notice 
that poles in regions III and IV have a real part extend- 
ing outside the band edges (indicated with an horizontal 
dotted line). However, as they enter region IV they loose 
their imaginary part. At this point the poles jump up to 
the second Riemann sheet, as represented in Fig. |4l and 
they become virtual states [i^l . This is manifested in the 
fact they are poles of Eq. [7] only because the unphysical 
branch of the self-energy was used. Hence, they do not 
show up as peaks or deltas in the density of states. In- 
deed, the LDOS given by Eq. [15] integrates to 1 within 
the band support. Only when the poles reach the band 
edge again, see Fig. [JJ they return to the first Riemann 
sheet to become localized states. 

Figure [5| displays the dynamical regions in the param- 
eter space. Notice that if we move along an horizontal 
line, which means keeping Vab constant, we are not going 
to see all the regions. At least one region would escape 
the analysis. From the experimental point of view this 
implies that, in order to find all the regimes by control- 
ling only one parameter, it is better to vary the system 
frequency than the system environment interaction. For 
that reason, in the plots that follow, we only vary the 
parameter Vab while keeping Vq constant at 0.8V^. 

In order to ensure consistency, we also studied the 
behavior of PAA{t) by an exact diagonalization of the 
Hamiltonian of the finite system. In this case the en- 
vironment size is taken large enough so the mesoscopic 
echoes do not show up at the times of interest [6lj. The 
evaluation of Paa (t) allows us to univocally identify the 
localized states and the different decay laws of the other 
regimes. In the next sub-sections we will show, in more 
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FIG. 4: (Color online) Paths in the complex energy plane of 
the two poles (black and red) of the GF as Vab decreases. 
Units are arbitrary. They go from the localized bonding and 
antibonding states into respective resonances that eventually 
collapse at the exceptional point. The bottom Riemann sheet 
contains the physical poles, while the upper sheet has the 
nonphysical poles. The broader horizontal line in the center 
represents the continuous band of environment states. Notice 
that localized states transform into virtual states and out-of- 
band resonances before becoming well defined resonances. 
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FIG. 5: Phase diagram of the different regimes as a function 
of the dimensionless hoppings Vab/V and Vo/V. The vertical 
dash-dotted line represents the condition used in Figs. (2] |4l 
E andU 



detail, the behavior of the system in each region. 



A. Region I: Collapsed resonances (overdamped 
decay) 



This region is found for: 



2 



\Vab\ < V-{V^-Vo^) 



(19) 



Under this condition, Paa (t) decays exponentially and 
without oscillations until the survival collapse time [6lj. 
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FIG. 6; (Color online) (a) Left side panels: (Black solid line) LDOS (in units 1/V) for the different regions [ Vab/V = 
0.35 (I), 1.0(11), 1.58 (III), 1.62 (IV), and 1.9 (V)]. (Blue dotted line) Li (e) and L2 (e). (Black dash-dotted line) Ni (e). The 
right side scale corresponds to the Li and L2 plots. (b) Right side panels: Survival probability Paa (t) as a function of time in 
h/V units. 



At this moment the return amplitude from the environ- 
ment starts to be comparable with the pure survival am- 
phtude interfering destructively [see Fig. [Stl-b)]. If we 
set Vp = IV we arrive to the case already treated in 
Ref. [6l| where it is analyzed the decay of a surface spin 
excitation when it interacts with a spin chain. A sim- 
ilar dynamics was recently predicted for a quasiparticle 
excitation at adsorbates in a surface band [62| . 

The real part of both poles of the GF coincide with 
the site energy, which means that the effective frequency 
is zero. However, their respective imaginary parts differ 
substantially (see Fig. [4]) . One of them moves away from 
the real axis as the S — £ interaction increases while the 
other approaches the real axis. This means that one state 



is captured by the environment while the other becomes 
isolated by cause of the quantum Zeno effect @. 

In this region Li and L2 are two LF centered at 0, but 
with different widths, 



[e) = C- 



2T 



1,2 



r2 ' 

1,2 



where 



C = 



1/2 T/2 

4rir2 ' 



(20) 



(21) 
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FIG. 7: Local dynamics in log-log representation for different 
regimes. Time is in YijV units. Curves II and V represent 
evolutions inside regions II and V, respectively. Curves II- 
III, III-IV, and IV-V represent the survival probability at the 
transitions between regimes. 



Figure [6l(I-a) shows the behavior of Na-, Ni and the 
Lorentzians functions: Li and L2. The centers and the 
linewidth of Li (e) and L2 (e) are exactly equal to the real 
and imaginary part of the poles of the GF, respectively. 

When Vab reaches zero, Fi vanishes, which is consis- 
tent with the fact that the first site becomes completely 
isolated and the second site behaves exactly as the case 
treated in Rcf. [6l|. On the other hand at the value 

2 _ T/2\l/2 



Va 



V- V' 



the system presents an EP. 



In this point both Lorentzians widths become equal. Be- 
yond this point lj starts to grow up (see Fig. [5]). We 
can interpret this critical point as a change of represen- 
tation from an adsorbed atom at surface of a semi- infinite 
chain, to an adsorbed dimer that appears when the initial 
adatom binds strongly to the surface atom. 



B. Region II: Resolved resonances (damped 
oscillations) 



This regime occurs for values of Vab satisfying 



V- V'-Vr 



2\l/2 



< \Vab \ < 




r 



(23) 



This parametric region is characterized by the oscillatory- 
exponential decay of Paa (t) at short times, Figs. [Ull-b). 
After the survival collapse it follows a decay controlled 
by quantum diffusion {t~^). This is better appreciated in 
Fig. [71 Indeed this regime was recently found in a model 
that describes the survival probability of an excitation 
in a multibarrier superlattice [63*1 . The power law decay 
is consequence of the environment's memory effects and 
corresponds to quantum diffusion from the "edge" of a ID 
system. This is verified by making a Fourier Transform 
of a LDOS of the form iVi(e) oc [e] which leads to 
a survival probability with an asymptotic form P(t) ^ 
^-(2i/+2)_ a bulk of d dimensions v = doff./2 — 1 with 
doff = d. Notably, our case corresponds to a survival 
probability of a state in a surface or edge, which decays 
as a bulk excitation in a higher effective dimension doff — 
d + 2n with n < d the order of the surface. In our case 
n = d = I. 



The poles in this region have both real and imaginary 
parts. At shorts times, the real part controls the oscil- 
latory behavior [lu — 2Re (e)] and the imaginary part 
determines the rate of the exponential decay. On the 
other hand, at long times the excitation decays with a 
t^'^ law and oscillates with a frequency determined by 
the environment band width 4]/. 

Here, the critical time tc is the time scale at which 
the quantum pathways returning from the environment 
starts to be comparable to the pure survival amplitude. 
Hence, tc, which decreases with Vab, divides the expo- 
nential decay from the quantum diffusive decay. 

The LDOS can be expressed in the same way as Eq. 
[T51 but in this regime the LFs are centered at symmetric 
points and have the same width [see Fig. (Sfll-a)], 



ii.2(e) = C 



2r 



(24) 
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with 



1/2 T/2 

ra _ ^0 ^ AB 



(25) 



It is interesting to notice that the Rabi frequency could 
either be slower than the unperturbed one, when Vab ^ 
V, or faster, when Vab ^ V. The precise cross over results 
from Eql23 when Er = Vab ■ 



Vo'-Vi^ {2V^-Vo^) 



4(^2 



WIb 

Hv^-v,')- 



vis {2V^-V,^)-Vi 
2{V^-Vi) 



(26) 



(27) 



C. Region III: Out-of-band resonances 
(environment controlled quantum diffusion) 

This region is found for Vab satisfying: 



'i + i6-^(y2_i/2) 



{2V^-VS)-2^V^ {V^^Vi)^ 



< \Vab \ < 



V- 



iv^-v^ 



1/2 



r 



(28) 
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FIG. 8: (Color online) Solid line: frequency u) evaluated from 
the poles distance (in units of V/h) vs Vab (in units of V). 
Blue dots: frequency fitted from the dynamics. The value for 
Vb is the same as in Fig. [2] 



If we only take into account the poles of the GF the 
system resembles the Resolved Resonances regime. How- 
ever, if we compare their LDOS and Paa (t), we conclude 
that this is a new regime. 

In this region, it is no longer possible to distinguish an 
exponential decay in Paa (t) [see Fig. [Hllll-b)]. This is 
because the time tc is too short. This implies that the 
oscillations are a direct consequence of the environment. 
In Fig. [8] it is observed that the frequency obtained di- 
rectly from the numerical solutions remains constant at 
4:V while uj (frequency evaluated from the poles) contin- 
ues to grow, confirming our analysis. 

The expression for the LDOS presented in Eq. [15] is 
still valid, but the LF are centered outside the band edges 



[see Fig. |ni[III-a)]. If we analyze the dynamics in term 
of the LDOS, we see that the tails of the LF decay as 
which means that the LDOS will have a Van Hove 
singularity as e^/^ x £~^, which implies that the local ex- 
citation will decay as at long times [see Figs. ini[IITb) 
andU]. 



D. Region IV: Virtual states (anomalous diffusion) 

This region corresponds to the range 
V+{V^-Vo^f^\ < \Vab\ < ^2{2V^-VS). (29) 

The poles in this region do not have imaginary part. For 
this reason it might be expected that poles were localized 
states dll. However, both the LDOS of Fig. E^IV-a) and 
the dynamics shown in Fig. [6lJIV-b) and Fig. [7] probes 
that this is not the case. The reason is that the poles 
had moved to a second Ricmann sheet [3, (see Fig. 
S]). In such case only the use of the unphysical sign of 
the self-energy could provide poles. The dynamics in 
this regime presents a striking transition between a 

decay [at Vab = V + {V'^ - Vi^f^ ] to a i"! behavior 



[at Vab = ^2 {2V^ - V§)\. This is shown in Fig. H 

Once again, we can express the LDOS as in Eq. [TSl 
however now the interpretation is different. In this region 
Sr is still outside the band edges and yet the value of F be- 
comes imaginary. The fact that F is transformed into an 
imaginary number implies that L\ and L2 are no longer 
LFs. If we now analyze the tails of these functions, we ob- 
serve that they decay as . This leads to a LDOS with 
a Van Hove singularity of the form e^^"^ x e^^. Therefore, 
we can achieve a behavior of Paa at long times. This 
fact is indeed confirmed by the observed dynamics (see 
Fig. [7]). As in the previous subsection, uj do not follows 
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the observed oscillation frequency [see Fig. [8] and [BIJIV- 
a)]. Instead, it is fixed by the environment bandwidtli 
AV. 

When Vab reaches the value \/2 {2V'^ — Vq), there is 
a change in the nature of the van Hove singularities from 
£^/^ to e~^/^. Consequently, PAA{t) decays exactly as a 
t~^{see Fig. [7]). From this point on, the states become 
localized. 

It is interesting to note that while the presence of the 
virtual states is not clearly distinguishable in the observ- 
able LDOS, the anomalous diffusion, where PAA{t) moves 
gradually between t~'^and t^^, should enable its experi- 
mental identification. 



E. Region V: Pure point states (localized) 

Finally, the last region appears when, 

\Vab\ > ^2 (2^2 -^0^). (30) 

Two localized states emerge from the band edges as 
shown in Fig. [6^V-a). The poles are real. Figure [8] 
shows that uj recovers its interpretation as the effective 
system frequency. In this region the environment renor- 
malization is almost negligible and its only effect is to 
slightly correct the value of the effective frequency. If 
Vab becomes large enough lu reaches ujq. The dynamics 
in this region. Fig. [Sl^V-b) is characterized by an oscil- 
latory PAA{t) that only decays for a short period after 
which the amplitude of the oscillation remains constant. 



V. CONCLUDING REMARKS 

By considering an exactly solvable model for a swap 
gate in presence of an environment, we discussed how 
the memory of the bath affects the dynamics beyond 
the Fermi golden rule. The unperturbed Rabi frequency 
sweeps through different dynamical and analytic regimes. 
Depending on the value of the internal time scales of the 
system, the environment can induce multiple phase tran- 
sitions in the system dynamics. Our model shows uncom- 
mon regimes as the exchange narrowing starting at the 
exceptional point and the virtual states by moving a sin- 
gle parameter. The fact that all the dynamical phases 
appear in the same system, offered the opportunity to 
study the transitions between them and hence to deter- 
mine the precise critical points. 

By studying the dynamics, we have characterized the 
difference between collapsed resonances, resolved reso- 
nances, out-of-band resonances, virtual states and the 
pure point states. In particular, we have shown that 
in the virtual state regime there is an anomalous quan- 
tum diffusive law, which at long times is observed as 
a change in the decay law from to t~^. We have 
found an expression for the LDOS where the presence 



of the resonances becomes explicit. This LDOS is fac- 
torized in three terms: a density of directly connected 
states (i.e. the LDOS at the environment's surface) and 
two Lorentzian Functions whose widths become imagi- 
nary in the virtual state regime. The edges of the LDOS 
determine the behavior at long times. It is then clear 
that the anomalous diffusion is related to the fact that 
the van Hove singularities at these edges are modified. 
This change occurs when the Lorentzian widths become 
imaginary. 

Note that the complexity of the dynamics for this 
relatively simple system, emerges as a consequence of 
the explicit way in which the environment is modeled. 
Details such as out-of-band resonant states and virtual 
states could not have been observed in simpler represen- 
tations of the environment as the usual broad-band or 
the self-consistent Born approximations. The results of 
our model system shows that a zero imaginary part of 
the poles is not enough as a localization criterion [58.] . In 
particular, virtual states have zero imaginary part, but 
any local excitation in this parametric regime shows a 
complete decay to the environment. 

Another parameter that characterizes the dynamics, is 
the oscillation frequency. First, in region I (collapsed res- 
onances), before the system reaches the EP, there is an 
overdamped decay as expected. When the unperturbed 
Rabi frequency ujq exceeds the critical value, at the EP, 
the resonances become resolved (region II) . Paa {t) show 
an exponentially attenuated oscillation with a frequency 
LU given by the poles difference. This frequency (D goes 
from 0, at the EP, to a value higher than ojq. In the re- 
gion III (out-of-band resonances) and IV (virtual states) 
the observed frequency is fully determined by the envi- 
ronment bandwidth and not by the GF poles, while the 
decay follows different power laws. These are particu- 
larly stable regimes for a swap gate since the effective 
Rabi frequency does not depend on the internal parame- 
ter. Further increase of leads to Region V of localized 
states where the observed frequency is always higher than 
ujo but tends to it as Vab — > oo. 

An interesting message that can be extracted from the 
phase diagram of Fig. [5] is that by changing Vab /V 
one can always move between all five regimes. However, 
for small values of the interaction with the environment 
(Vq/V) the transition to localized states occurs within a 
very small range of Vab /V and hence most chances are 
that it goes unobserved. The diagram indicates that a 
clear numerical or experimental observation of the full 
dynamical wealth would only be possible for Vq/V < 1. 

One might wonder how decoherence affects the swap 
gate. This, of course, would depend on the specific situa- 
tion. Consider, for example, a spin chain with XY inter- 
action as described in Refs. [l^l and [9]. One may assume 
that the initial condition of the dimer used a swap gate 
is in the state |1,0) while the chain that acts as environ- 
ment is in a mixed state with 1/2 occupation at each site. 
The environment density will gradually tunnel into the 
dimer until it becomes a mixed state itself. In this case. 
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decoherence originates in the mixed initial condition of 
the "bath" spins. While the probability would tends to 
an asymptotic value, which might differ from the ergodic 
1/2 64] , the Rabi frequency would remain that evaluated 
from our coherent calculation. Another situation that 
one might conceive is a second environment £' yielding 
a weak decoherent process of strength 77 added homoge- 
neously throughout the whole system as an imaginary 
component of the site energies. In this case, oscillations 
would simply degrade at a rate 2ri/h. However, a pre- 
diction of the effect of strong decoherence requires spe- 
cific analysis as it is known that, even at a FGR level, 
self-consistent decoherence may lead to dynamical phase 
transitions [13, Hi]. 

As a final remark, we again mention that because of its 
simplicity, our model can be arranged to describe various 
physical systems beyond spin chains such as microwave 
antenna arrays [gFI, energy shuttling on plasmon wires 
made of nanoparticles 66], arrays of tunneling coupled 
optical waveguides [67|, periodic elastic arrays [gl], or 
acoustic time reversal cavities [g^ . Reciprocally, most of 



them should present the variate dynamical phenomena 
discussed here, provided that one focuses in the proper 
parameter range. These examples also suggest possible 
experimental setups where the parameters found in this 
work can be used as a knob enabling to store and ex- 
change energy. 



Acknowledgments 

The authors acknowledge Gonzalo Alvarez, Hernan 
Calvo, Ernesto Danieli, and Elena Rufeil-Fiori for shar- 
ing their experiences in related fields. Patricia Levstein 
and Pablo Serra helped with discussions and useful refer- 
ences. HMP acknowledges hospitality of MPI-PKS and 
Abdus Salam ICTP where discussion with Ingrid Rotter 
and Maksim Miski-Oglu helped to enrich this work. The 
work was made possible through the financial support 
from CONICET, ANPCyT, and SeCyT-UNC. 



[1] Charles H. Bennett and David P. DiVincenzo, Nature 

(London) 404, 247 (2000). 
[2] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 

K. C. Nowack, T. Meunier, L. P. Kouwenhoven and L. 

M. K. Vandersypen, Nature (London) 442, 766 (2006). | 
[3] D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. 

Urbina, D. Esteve and M. H. Devoret, Science 296, 5569 

(2002). 

[4] Y. Nakamura, Yu. A. Pashkin and J. S. Tsai, Nature 

(London) 398, 786 (1999). 
[5] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. 

Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. 

Schoelkopf, Nature (London) 431, 162 (2004). 
[6] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 

Phys. 73, 565 (2001). 
[7] F. Schmidt-Kaler, H. Haffner, M. Riebe, S. Guide, G. 

P. T. Lancaster, T. Deuschle, C. Becher, C. F. Rocs, J. 

Eschner and R. Blatt, Nature (London) 422, 408 (2003). 
[8] C. Ramanathan, N. Boulant, Z. Chen, D. G. Cory, C. 

Isaac and M. Steffen, Quantum Inf. Processing 3, 1 

(2004). 

[9] G. A. Alvarez, E. P. Danieh, P. R. Levstein, and H. M. 

Pastawski, J. Chem. Phys. 124, 1 (2006). 
[10] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003); W. H. 

Zurek, F. M. Cucchietti, and J. P. Paz, Acta Phys, Pol, 

B 38, 1685 (2007). 
[11] C. J. Myatt, B. E. King, Q. A. Turchette, C. A. Sack- 

ett, D. Kielpinski, W. M. Itano, C. Monroe and D. J. 

Wineland, Nature (London) 403, 269 (2000). 
[12] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 

756 (1977). 

[13] S. Pascazio and M. Namiki, Phys. Rev. A 50 4582 (1994). 
[14] H. M. Pastawski and G. Usaj, Phys. Rev. B 57, 5017 
(1998). 

[15] J. H. Van Vleck, Phys. Rev. 74, 1168 (1948). 

[16] N. Bloembergen, E. M. Purcell, and R. V. Pound, Phys. 



Rev. 73, 679 (1948). 
[17] P. W. Anderson, J. Phys. Soc. Jpn. 9, 316 (1954). 
[18] G. A. Alvarez, E. P. Danieh, P. R. Levstein, and H. M. 

Pastawski, Phys. Rev. A 75, 062116 (2007). 
[19] E. P. Danieh, H. M. Pastawski, and G. A. Alvarez, Chem. 

Phys. Lett. 402, 88 (2005). 
[20] J. R. Taylor, Scattering Theory: The Quantum Theory 

of Non-Relativistic Collisions (Dover, Mineola, 2000), p. 

246. 

[21] H. M. Nussenzveig, Causality and Dispersion Relations 
(Academic, New York, 1972) pp. 158,159, 220. 

[22] M. K. Oberthaler, R. Abfaherer, S. Bernet, J. Schmied- 
mayer, and A. Zeilinger, Phys. Rev. Lett. 77, 4980 
(1996). 

[23] T. Stehmann, W. D. Heiss, and F. G. Scholtz, J. Phys. 

A 37, 7813 (2004). 
[24] A. L. Shuvalov and N. H. Scott, Acta Mech. 140, 1 

(2000). 

[25] M. V. Berry and M. R. Dennis, Proc. R. Soc. London, 

Ser. A 459, 1261 (2003). 
[26] G. Khitrova, H. M. Gibbs, M. Kira, S. W. Koch, and A. 

Scherer, Nat. Phys. 2, 81 (2006). 
[27] C. Dembowski, H. D. Graf, H. L. Harney, A. Heine, W. 

D. Heiss, H. Rehfeld, and A. Richter, Phys. Rev. Lett. 

86, 787 (2001). 

[28] B. Dietz, T. Friedrich, J. Metz, M. Miski-Oglu, A. 

Richter, F. Schafer, and C. A. Stafford, Phys. Rev. E 

75, 027201 (2007). 
[29] C. Dembowski,B. Dietz, H. D. Graf, H. L. Harney, A. 

Heine, W. D. Heiss, and A. Richter, Phys. Rev. E 69, 

056216 (2004). 

[30] C. Dembowski, B. Dietz, H. D. Graf, H. L. Harney, A. 

Heine, W. D. Heiss, and A. Richter, Phys. Rev. Lett. 90, 

034101 (2003). 
[31] G. Weinreich, J. Acoust. Soc. Am. 62, 1474 (1977). 
[32] L. M. B. Napolitano, O. R. Nascimento, S. Cabaleiro, J. 



11 



Castro, and R. Calvo, Phys. Rev. B 77, 214423 (2008); A. 

J. Costa-Filho, C. E. Munte, C. Barberato, E. E. Castel- 

lano, M. P. D. Mattioli, R. Calvo, and O. R. Nascimento, 

Inorg. Chem. 38, 4413 (1999); Rafael Calvo, Appl. Magn. 

Reson. 31, 271 (2007), and references therein. 
[33] C. A. Stafford and B. R. Barrett, Phys. Rev. C 60, 

051305(R) (1999). 
[34] W. D. Heiss and A. L. Sannino, Phys. Rev. A 43, 4159 

(1991). 

[35] F. Keck, H. J. Korsch, and S. Mossmann, J. Phys. A 36, 

2125 (2003). 
[36] W. D. Heiss, Phys. Rev. E 61, 929 (2000). 
[37] M. Miiller, F. M. Dittes, W. Iskra and I. Rotter, Phys. 

Rev. E 52, 5961 (1995). 
[38] F. Stefani and G. Gerbeth, Phys. Rev. Lett. 94, 184506 

(2005) . 

[39] D. M. Cardamone, C. A. Stafford, and B. R. Barrett, 

Phys. Stat. Sol. B 230, 419 (2002). 
[40] E. P. DanieU, G.A. Alvarez, P. R. Levstein, and H. M. 

Pastawski, Solid State Commun. 141, 422 (2007). 
[41] S. Longhi, Phys. Rev. A 74, 063826 (2006). 
[42] H. L. Calvo and H. M. Pastawski, Braz. J. Phys. 36, 963 

(2006) . 

[43] H. A. Weidenmiiller, Phys. Rev. B 68, 125326 (2003). 

[44] H. Hogreve, Phys. Lett. A 201, 111 (1995). 

[45] B. Marcelis, E. G. M. van Kempen, B. J. Verhaar, and S. 

J. J. M. F. Kokkelmans, Phys. Rev. A 70, 012701 (2004). 
[46] A. M. Pupasov, B. F. Samsonov, and J.-M. Sparenberg, 

Phys. Rev. A 77, 012724 (2008). 
[47] P. Serra, S. Kais and N. Moiseyev, Phys. Rev. A 64, 

062502 (2001). 

[48] M. T. Yamashita, T. Frederico, A. Delfino and Lauro 

Tomio, Phys. Rev. A 66, 052702 (2002). 
[49] Horacio M. Pastawski and Ernesto Medina, Rev. Mex. 

Fis. 47, 1 (2001) cond-mat/0103219 
[50] E. Danieli, H. Pastawski and P. Levstein, Chem. Phys. 

Lett. 348, 306 (2004). 
[51] Z. L. Madi, B. Brutscher, T. Schulte-Herbriiggen, R. 

Briischweiler, and R. R. Ernst, Chem. Phys. Lett. 268, 

300 (1997). 



[52] S. Chakravarty and A. J. Leggett, Phys. Rev. Lett. 52, 
5 (1984). 

[53] M. C. Gutzwiller, in Chaos and Quantum Physics, edited 
by M.-J. Giannoni, A. Voros, and J. Zinn- Justin (North- 
Holland, Amsterdam, 1991). 

[54] Paolo Facchi and Saverio Pascazio, La Regola d'Oro di 
Fermi (Bibliopolis, Napoh, 1999). 

[55] L. Khalfin, Sov. Phys. JETP 6, 1053 (1958). 

[56] L. Fonda, G. C. Ghirardi, and A. Rimini, Rep. Prog. 
Phys. 41, 588 (1978). 

[57] G. Garcia-Calderon, J. L. Mateos, and M. Moshinsky, 
Phys. Rev. Lett. 74, 337 (1995). 

[58] E. N. Economou, Green's Functions in Quantum Physics, 
3rd ed. (Springer, Heidelberg, 2006). 

[59] H. M. Pastawski, Physica B 398, 278 (2007). 

[60] E. Narevicius, P. Serra, and N. Moiseyev, Europhys. Lett. 
62, 789 (2003). 

[61] E. Rufeil-Fiori and H. M. Pastawski, Chem. Phys. Lett. 

420, 35 (2006). 
[62] B. Gumhalter, A. Siber, H. Buljan, and T. Fauster, Phys. 

Rev. B 78, 155410 (2008). 
[63] G. Garcia-Calderon, R. Romo, and J. Villavicencio, 

Phys. Rev. B 76, 035340 (2007). 
[64] H. M. Pastawski, P. R. Levstein, and G. Usaj, Phys. Rev. 

Lett. 75, 4310 (1995). 
[65] Luke A. Sweatlock, Stefan A. Maier and Harry A. At- 

water. Proceedings of the 53rd Electronic Components 

and Technology Conference, New Orleans, 2003, (IEEE, 

Washington D.C., 2003), p. 1648. 
[66] M. L. Brongersma, J. W. Hartman, and H. A. Atwater, 

Phys. Rev. B 62, R16356 (2000); R. Baer, K. Lopata and 

D. Neuhauser, J. Chem. Phys. 126, 014705 (2007). 
[67] S. Longhi, Phys. Rev. Lett. 97, 110402 (2006). 
[68] L. Gutierrez, A. Dfaz-de-Anda, J. Flores, R. A. Mendez- 

Sanchez, G. Monsivais, and A. Morales, Phys. Rev. Lett. 

97, 114301 (2006). 
[69] H. L. Calvo, E. P. Danieh, and H. M. Pastawski, Physica 

B 398, 317 (2007). 



